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1 Introduction 



Let K s be the trace-class operator with kernel 

Ai(t)Ai'(tQ-Ai( M )Ai'(t) 
K s (t,u) = ^— — (1) 

(see [3 1 J ) acting on L 2 (—s, oo). Here Ai (x) is the Airy function (see, e.g., pQ). In this paper 
we are concerned with the behavior of det(J — K s ) as s -> +oo. Our main result is the 
following. 

Theorem 1. The large-s asymptotic behavior of the Fredholm determinant det(J — K s ) is 
given by the formula 

lndet(I - K.) = ~ - Uns + X + 0{s^ 2 ), (2) 

where 

X = ^ln2 + C / (-l), (3) 

and is the Riemann zeta- function. 

The Airy-kernel determinant det(J — K s ) is the edge scaling limit for the largest eigen- 
value of a random n x n Hermitian matrix H from the Gaussian Unitary Ensemble (GUE) 
(see [291 E]) as n -> oo: More precisely, if \\(H) > X 2 (H) > • • • > A n (,ff) denote the 
eigenvalues of H, then 

det(J - K s ) = lim Prob {iJ G GUE : (Ai(#) - v / 2^)2 1/2 n 1/6 < -s} (4) 

n— >oo 

(See pT] . [31], and also [8] for some history of (J3J)). 

This determinant also describes the distribution of the longest increasing subsequence 
of random permutations [3l [25]. Namely, let tt = ix%% ■ ■ - in be a permutation in the group 
S n of permutations of 1, 2, . . . , n. Then a subsequence ik 1 ,h 2 i ■ ■ -h r i hi < ^2 < • • • < K, of 
7r is called an increasing subsequence of length r if < i& 2 < • • - < Zfc r . Let Z„(7r) denote 
the length of a longest increasing subsequence of 7r and let S n have the uniform probability 
distribution. Then l n (ir) is a random variable, and [3] 

det(J - K s ) = lim Prob {n e S n : (l n (n) - 2^l)n- 1/6 < -s} (5) 

The distribution Ftw{x) = det(J — K_ x ), known as the Tracy- Widom distribution, 
admits the following integral representation [31] : 

(y-x)u 2 (y)dy \, (6) 



F TW (x) = exp 
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where u(y) is the (global) Hastings-McLeod solution of the Painleve II equation 

u"(y)=yu(y) + 2u 3 (y), (7) 
specified by the following asymptotic condition: 

u(y) ~ Ai(y) as y — > +00. (8) 
The behavior of u(y) as y —>■ —00 is given by the relation [22] : 



n(y) = J-|(l + -i I + 0(y- 6 ) ), y^-oc, (9) 



2 V % 3 
from which one learns that as s — > +00, 

r°° s 3 1 

-J (s + y)u 2 (y)dy + — + -\ns = as + b + o(l) (10) 

for some constants a, b. The content of Theorem 1 is that a = and b = x as in ([3]). The 
value ([3]) of the constant x was conjectured by Tracy and Widom in [31] on the basis of the 
numerical evaluation of the l.h.s. of (1101) as s — > +00 and by taking into account the Dyson 
formula for a similar constant in the asymptotics of the so-called sine-kernel determinant 
[3T] . The sine-kernel determinant describes the gap probability for GUE in the bulk scaling 
limit as n — > 00 [29|. 

Dyson's conjecture for the constant in the asymptotics of the sine-kernel determinant 
was proved rigorously in independent work by Ehrhardt [TH] and one of the authors [26J, and 
a third proof was given later in [T7j. The two latter works use a Riemann-Hilbert-problem 
approach. The proof in [26] relies on a priori information from [33], whereas the proof in 
[T7] is self-contained. The proof of Theorem 1 in this paper follows the method in [17| . 

As discussed in [17], the key difficulty in evaluating constants such as x in |2} in the 
asymptotic expansion of the determinants, is that in the course of the analysis one most 
naturally obtains expressions only for the logarithmic derivative with respect to some aux- 
iliary parameter, say a, in the problem, and not the determinant itself. After evaluation of 
these expressions asymptotically, the constant of integration remains undetermined. In [17] 
and [26J, this difficulty is overcome by utilizing a scaling limit of finite-ra random matrices 
together with universality in the sense of random matrix theory (see, e.g., [15]), in a way 
that is inspired by, but different from, Dyson [IB]. We proceed as follows. 

Consider the scaled Laguerre polynomials Pk{x) defined for some integer n by the or- 
thogonality relation 



e nx p k (x)p m (x)dx = 6 ktTn , fc,m = 0,l..., (11) 

The polynomial pk(x) = xux k + • • • is of degree k and is related to the standard Laguerre 
polynomial Lf\x) (see [30]) as follows: 



p k (x) = 2v^4 0) (4 



nx 
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with leading coefficient 



^ = (-l)^(4n) fc . (12) 

The scaling here is chosen so that the asymptotic density of zeros of the polynomial p n (x) 
(with index n) is supported on the interval (0,1) (as opposed to (0, An) for Ln(x)). See 
E01 M and below. 



In the unitary random matrix ensemble defined by the Laguerre weight, the distribution 
function of the eigenvalues is given by the expression: 

- n— 1 

dP{x , . . . , x„_ x ) = — — I j (x t - Xjf Y[ e- 4x > n dxj, (13) 

0<i<j<n-l j=0 

where the normalization constant 

i poo poo n ~~ 1 

c » = ^ / ■■■ n - n ' '" v/ -'> ( i4 ) 

' ^° ^° 0<i<j<n-l j=0 

By a well known identity (e.g. [30], [7]), the r.h.s. of the above expression gives 

n—1 n—1 

^ = n^ 2 = (4n)-" 2 n^ (15) 

fc=0 fc=0 

where f lT2|) was used. 

For a > 0, the probability D n (a) that the interval (a, oo) has no eigenvalues is given by 



dm =-i—( a ... r n ( Xi - X] f n e -^ n d Xjl ^ 

""' ^° 0<i<i<n-l j=0 



o Jo 



By standard arguments (cf. [3[29]), this quantity can be written as the Fredholm determi- 
nant of an integral operator on £2(0, 00) in the following way: 

!>„(«) = det(J - K n X(a,oo)), ^n^i V) = T , ( 17 ) 

where 

uj k (x) = e~ 2nx p k (x) } fc = 0,l,..., (18) 
and X(a,oo) is the characteristic function of the interval (a, 00). 
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If x = 1 + l/(2n) + -u/(2n) 2 / 3 with it fixed, then as n — > oo, we obtain from classical 
results on the asymptotics of the Laguerre polynomials (see [32l 130] ) : 



Wn 1 + TT" + 



1/ 



2n (2n) 2 / 3 



u; n „i 1 + — + 



2n (2n) 2 / 3 
0(n- 2 / 3 )]) = (-l)'- 1 



(i-(4n + 2 + 2(2^1*)) = (-lf^ {Ai («) + 0(rT 2 / 3 )} ; 
4(n - 1) + 2 + 2(2[n - 1]) 1/3 ( it + 



^ 1 4^ 



(2n)V3 



2^ 



Ai U + 



(2n)V3 



+ 0(n 



-2/3>! 



(2[n- 1])V3 

where Ai (x) is the standard Airy function. Let 

Ai {u)Ai'{v) - Ai (u)Ai'(it) 



(19) 



■K a i r y (it, 



U — V 



Set 
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(n) 



1 + 



1t 



+ 



(2n) 2 / 3 2n 

It follows from (jl~9l) that for any fixed it, u we have 

1 



1 

+ — . 



(2n) 2 / 3 2n 



lim , 

n-oo (2n) 2 / 3 



K n { 



(«) An) 



airy V ) • 



In fact, this asymptotics is uniform for it, v > Lq, where Lo is an arbitrary constant. Indeed, 
for any L there exists C = C(L ) > 0, c = c(L ) > such that 



&>d k 



r^ n (u (n) ,i/ n) ) -# oin/ (it,i;) 



< C- 



n 



2/3 



(20) 



* w V(2n) 2 / 3J 

u,v>L , j, k = 0, 1. 

This estimate can be proved in a same manner as estimate (3.8) in [8]. In [8] the authors use 
global estimates for orthogonal polynomials on M taken from [IB] : Here the relevant global 
estimates can be obtained from [32] . 

As in [Sj, estimate (120]) immediately implies that for any fixed sGl, 



lim D 7 

n^oo 



(21) 



Below we obtain the asymptotics of the determinant det (/ — K aiTy X(- s ,oo)) = det (J — K s 
as s — > +oo. In order to do this, we analyze the asymptotics of (1T7I) for all a from a close 
to zero to a = 1 — s/(2n) 2 / 3 . Note that the determinant (1T7I) has the structure of so-called 
integrable determinants [23]. Therefore, it is not surprising that there exists a differential 
identity for lnD n (a) in terms of the solution of a related Riemann-Hilbert problem. Solv- 
ing the Riemann-Hilbert problem asymptotically as n — > oo, we find the asymptotics of this 
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logarithmic derivative uniform for a £ [0, 1 — so/(2n) 2,/3 ], (2n) 2 / 3 > so for some (large) so > 0. 
Integrating these asymptotics from a close to zero to a = 1 - s/{2n) 2/3 , s Q < s < (2n) 2 / 3 
we obtain the asymptotics of D n (l — s/ (2n) 2 ^ 3 ) provided we know the asymptotics of D n (a) 
for a close to zero. The latter, however, is readily obtained from the series expansion of the 
multiple integral formula for D n (a) (see (I22II27P below). More precisely, the "inner workings" 
of the method in this paper (cf. also (133) in [17]) can be seen from formula (11611) below, 
which is obtained by integrating the derivative (d / da') In D n (a') from a' = cto to a' = a. 
The key fact is that the estimate on the derivative is uniform for < a' < 1 — s/(2n) 2//3 , 
s > so (see (115211153!) ): This leads to the error estimate 0(l/(n(l — a) 3 / 2 )) in (11611) . Using 
( 1271) . we can then let «o — * 0: The singularities on the l.h.s. and the r.h.s. of (11611) cancel 
out, and we are left with (11621) . Using (11621) . we immediately obtain Theorem 1. Note that 
in our calculations formula ([3]) for \ does not arise from an evaluation of D n (ao) as n — > oo 
for some fixed «o- Rather it arises, somewhat paradoxically, from the behavior of D n (ao) as 
a — > with n fixed as given in (1271) . 

In Section 2 the series expansion of D n (a) for n fixed and a — > is derived, as indicated 
above. In Section 3 we obtain an asymptotic (n — > oo) solution of the Riemann-Hilbert 
problem related to (1T7|) . Moreover, in Section 3, a differential identity for ^ln.D n (a) is 
obtained in terms of the matrix elements (and their first derivatives) of the solution to the 
Riemann-Hilbert problem at the point a. An alternative derivation of this identity, which 
is closer to the spirit of integrable systems and r-functions (see, e.g., [3 El E3]), is given in 
the Appendix. The identity is then evaluated asymptotically in Section 4 using asymptotics 
found in Section 3. In Section 5 the identity is integrated, and the results of Section 2 are 
then used to complete the proof. 

Remark. Universality allows for considerable freedom in the choice of the approximating 
ensemble in the above method. We choose to consider the Laguerre ensemble, although we 
could have considered, for example, GUE itself: for GUE, however, the analysis turns out 
to be algebraically more complicated. (For example, in the GUE case there will be two 
endpoints instead of one endpoint at x — 1, see (fl2l) et seq.) In choosing the approximating 
ensemble, it is essential that the various constants that arise can be evaluated explicitly as 
in (12"7|) and also in formula (17) in [17J. In both cases we see that ultimately the formula for 
the desired constant arises from classical formulae for the Legendre polynomials. 

In physics, and also in mathematical physics, universality is often viewed as a passive 
statement that certain systems "behave in a similar fashion" . The thrust of this paper, going 
back to Dyson [18J, is that universality can be used as an active analytical tool to obtain 
estimates for asymptotic problems of mathematical and physical interest. 

Addendum. We draw the attention of the reader to the work [2] of Baik, Buckingham, and 
DiFranco, in which the authors give a different proof of ([3]) together with related results for 
GOE and GSE. The paper [2J appeared after our paper was written and refereed. 
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2 Expansion of D n (a) as a — ► 0. 

In this section we derive a series expansion for D n (a) as a — > 0. Changing the variables 
Xj = (a/2)(tj + 1) and expanding the exponent in ffTBT) . we obtain for fixed n: 



D n (a) = — (I ) / • ■ ■ / J] (f, - t,) 2 - 2«n(t, + 1) + 0(a 2 ))^. 

n?7 " ^ _1 0<?<j<n-l j=0 

1 / a \ " 2 

— (-) A n (l + 0„(a)), 

f22i 



where 

n-l 



a - = ^ / ••• / n (*i-*i) a n^ 

• J-l J-l <i<j< n _i 7= o 



(23) 



can be expressed in terms of the product of the leading coefficients (cf. (11411150 ) of the 
Legendre polynomials: 



^-ll[(2fc)!p2fc + l- (24j 

The asymptotics of A n as n — > oo (used first by Widom in [33], and then in [17]) are given 
by the expression 

In A n = -n 2 ln2 + nln(27r) - ^lnn+ — ln2 + 3C'(-1) + S n , n -> oo. (25) 

where C'( x ) is the derivative of Riemann's zeta-function, and 5 n — > as n — > oo. The 
zeta-function originates from the expansion of the product of factorials. 

The asymptotics of C n (|T5|) have a similar form, 



lnC n = -^- + ln4jn 2 + nln(27r)--lnn + 2C / (-l) + 5„, 5 n -> 0, n -> oo. (26) 

Substituting the asymptotics (I25II26P into ( |22l) . we obtain for a > 0: 

In D n (a) = Q + In a) n 2 - i In - + C'(-l) + 5 n + O n (a), (27) 

where 5 n depends on n only and 5 n — > as n — > oo. Note for later application (see proof of 
Lemma 2) that the error term O n (a) is analytic in a, in particular, (d/da)O n (a) = O n (l). 
We shall use formula (1271) in the last section. 



Caveat: O n (a) — > as a — > 0, n fixed: no claim is made about O n (a) as n — > oo. 
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3 Differential identity and the Riemann-Hilbert prob- 
lem 



3.1 Initial transformations 

In what follows, unless explicitly stated otherwise, we will always assume < a < 1. At 
certain points in the text, however, we will also consider a in a small neighborhood T> £0 (0) 
of a = (see the discussion in the end of Section 3.1.) 

The multiple integral (fl6|) can be written as (cf. (I14|ll5j) ): 

D n (a) = ^nV, ( 28 ) 

^n ■ n 
3=0 

where 6j are the leading coefficients of the polynomials qj(x) = 8jx 3 + • ■ ■ satisfying 

qk{x)q m (x)e~ 4nx dx = 5 km , k, m = 0, 1, . . . (29) 



It is convenient to write this orthogonality relation in the form 



qi {x)x k e- 4nx dx= 5 -^, fc = 0,l,...,j, j = 0,1,2,... (30) 

U 3 



Note, in particular, that 
d 



J 9i(^)^gj(^)e" 4na; rfx 



/ qAx) \^—i-x 3 + polynomial of degree less than j ) e 4nx dx — 
Jo \da J 



9j da 



Using relation ([3"Tj) . we obtain 



« , n— 1 n— 1 _, 7n n—1 

— In D n (a 
da 

3=0 3=0 J 3=0 



3=0 3=0 3 3=0 J U 

r a r> f n ~ l \ rl ( r aTl ~ x \ 71-1 



d=0 J V u 3=0 J 3=0 



(31) 



(32) 



By (1291) with k = m = j, the last integral (inside the brackets) in (|32|) equals n and 
hence vanishes upon differentiation. 

Applying the Christoffel-Darboux formula, 

n-1 ^ 

Yl = -^-Un^Mn-iix) - g n (x)g^_ 1 (x)), (33) 

3=0 n 



to the last sum in (1321) . we obtain 

±]nD n (a) = ^(^(aJ-^Ua))- (34) 
da u n 

Here and below the prime denotes differentiation w.r.t. the argument x. 

Formula (I3"l"|) shows that ^ln.D n (a) depends only on q n , q n -\. This property is crucial 
for the analysis below. 

As noted in [20], orthogonal polynomials can be represented in terms of a solution to 
an associated Riemann-Hilbert problem. In the present case, the relevant Riemann-Hilbert 
problem is formulated as follows: Find a 2 x 2 matrix-valued function V(z) satisfying the 
conditions: 



(a) V(z) is analytic for z£C \ [0, a}. 

(b) Let x G (0, a). V(z) has L 2 boundary values V+{x) as z approaches x from above, and 
V-(x), from below. They are related by the jump condition 

/ i -Anx \ 

V + (x) = V.(x)l Q \ ), xe(0,a). (35) 

(c) V(z) has the following asymptotic behavior as z — > oo: 

V(z) =(^ + o(j^J z™\ where a 3 = ( J . (36) 

This Riemann-Hilbert problem (RHP) has a unique solution for any n, a > 0, and, in 
particular, Vu(z) = q n {z)/6 n and V>zi{z) = —2i\id n _xQn-\{z)- Therefore we can rewrite the 
differential identity (!34j) in terms of Viz) in the formo 

^-\nD n (a) = ^l(Vu{a)V^{a) - K(a)V 21 {a)) . (37) 
da 2m 

In this section our task is to solve the RHP for V(z) asymptotically (in other words, 
to find asymptotics of the polynomials qk(z)) as n — >• oo. The results will then be used in 
Section 4.2. to evaluate the r.h.s. of (1371) . 

Following the steepest descent method for RH problems as described in [TBI E] , we first 
of all need to find a so-called ^-function: In the present situation this reduces to finding a 
function analytic outside the interval (— oo, a) and continuous up to the boundary with the 
properties: 

(a) g(z) = ln(z) + 0(1/ z) as z — » oo; 



1 An alternative derivation of this identity is presented in the Appendix. 
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(b) there exists a constant I such that the boundary values g±(x) = lim £ j g{x ±ze) of g(z) 
are related as follows: 

g+(x) + g-{x) - Ax - I = 0, xe(0,a); (38) 

(c) on (0, a), (7+(s) — gr_ (rc) is purely imaginary, and i(d/dx)(g + (x) — g~(x)) > 0; 

(d) e s+(*)-s-W = 1 on (-oo,0). 

A standard computation shows that if such a function g(z) exists than it is uniqueU 
Formally, the derivative g'(z) of g(z) must have the properties: 

(a') g'(z) = 1/z + 0(l/z 2 ) as z -> oo; 
(b') £f+(x) + (/L(a:) = 4 for x G (0, a). 

It is easy to verify that the following function satisfies these conditions: 

... 1 + a — 2z , . 

g '(z) = 2 + ^==. (39) 

(In fact, g'(z) is the unique function with L p boundary values g'± satisfying (a') and (b') 
for any 1 < p < 2.) In fl39l) . the branch is chosen so that a/ z(z — a) is analytic in the 
complement of (0, a) and positive for z > a. 

Therefore, 



g(z)= [ Z g'(t)dt + C, 

J a 



where the constant C is determined from the condition that g(z) — ln(z) = 0(l/z) as z — > oo. 
This gives 

g(z) = 2z - a + In - + / dt, (40) 



4 



a v/t(T 



and it is easy to verify that g(z) indeed satisfies (a)-(d). From (13811401) we now see that 

Ol 

/ = -2a + 21n-. (41) 
4 

We need to analyze the RHP for V(z) asymptotically as n — > oo uniformly for < a < 
1 — so/(2n) 2//3 where so is a fixed (large) number. The steepest descent method continues 
with the following steps (see [T6| [7]): 

1) the RHP for V is conjugated by e n9iz)a3 ; 



2 Note that as the contour for the RHP is (0,a), the extra condition (4.14) for g(z) in [16] is redundant 
in the present situation. 
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" 1 



Figure 1: Conformal mapping. 



2) the contour (0, a) is split into lenses; 

3) matching parametrices for the solution to the RHP are constructed (i) away from the 
end-points and a, (ii) in neighborhoods of and a, respectively. 

By means of these steps, the RHP reduces as n — > oo to a small norm problem which 
can be solved by a Neuman series. 

All these steps go through in the standard way except for the construction of the 
parametrix in a neighborhood of a. As we see from [T6| [7j the method requires that in 
a neighborhood \z — a\ < e, e small and fixed, 

(g + - gj){z) = (z- af{c + 0{z - a)), (42) 

for some c ^ and some exponent (3 > 0. (In [T5], (3 = 3/2.) In our case for < z < a, 

(g+ - g-){z) = 2 f 1 t^Z^l dt = -i=(2 _ a )V2(i _ Q + (z - a)). (43) 
J a yjt{t - a) V" 

For any fixed < a < 1 we see that (g + — g-)(z) satisfies (]4"2l) . As a — > 1, we have to 
make the neighborhood |z — a\ < e smaller and smaller. The constant c in (1421) depends 
then on a, but that, in itself, is not an insurmountable problem. The real problem is that, 
unlike the situation in [26J, the parametrix away from the points 0, a (see [I~6l [7]) contains 
certain terms of the form (z/(z — «)) 1//4 evaluated on {z : \z — a\ = e}, and as a result is 
not uniformly bounded when 1 — a, and hence e, approach zero. At the same time, there is 
not enough decay in the other relevant quantities to compensate for this. The problem can 
be circumvented, however, by introducing a transformation of the z-plane that "regularizes" 
the RHP in a neighborhood of z — a. Namely, set 

A= 1_ V ' z ^ 1. (44) 
a 1 — z 
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This fractional-linear transformation maps the interval [0, a] onto [0,1], the point z — 1 is 
mapped to infinity, and infinity is mapped to A = —(1 — a)/ a. The inverse transform is 



aX 

1 — a + aX ' 



1 — a 



a 



(45) 



Thus z(X) is analytic from C \ { — (1 — a)/a] into C, taking the complement of [0, 1] onto 
C\[0,a]. 

The fact that in our case we could not obtain an estimate of the form (1421) uniformly 
as a | 1 originates in the vanishing of the numerator in the integral for g + — g_ in (1431) at 
the point t = (1 + a)/2 G (a, 1). Under the transformation z — > A the point (1 + ct)/2 is 
mapped to A = 1 + or 1 . This point is at a positive distance from the contour < A < 1 for 
a G (0, 1). This means that we will be able to construct a parametrix for the solution of the 
RHP in the A variable in a fixed neighborhood about A = 1. On the other hand, the point 
A = —(1 — a)/a (the image of ^-infinity) now approaches the contour as a f 1, and we will 
need to contract the neighborhood of A = so that this point remains outside. We shall 
see, however, that this neighborhood presents no problem, as the relevant terms of the jump 
matrix for the final i?-RHP (see ( ITU]) and the argument after (11031) below) decay sufficiently 
fast on the boundary of the neighborhood. 

For any A G C\ ([0,1] U {-(!-«)/«}) set 



17(A) = V(z(\)), 

where z(X) = aX/ (1 — a + aX) as in fj4"5l) . 

Then we obtain the following Riemann-Hilbert problem for ^(A): 



(46) 



(a) U(X) is analytic for A G C \ ([0, 1] U {-(1 - a) /a}). 

(b) Let A G (0, 1). U has L 2 boundary values U + (X) as A approaches the real axis from 
above, and t7-(A), from below. They are related by the jump condition 



U + {X) = U-(X) 



e -4nz(A) 
1 



AG (0,1). 



(47) 



(c) U(X) has the following asymptotic behavior as A — » — (z — > oo): 

1 



U(X) 



7 + 



z(X) 



z(X) 



(48) 



We transfer g{z) to the A-plane by defining 



g(X)=g(z(X)), for AgC\ 



1 — a 



a 



(49) 
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Necessarily, g(X) is analytic on its domain. We obtain 

-/an n r i\ 1 — ^ i a \s/2 /" A 1 + Qt(l— t) (it 

S"(A) = 2a (a - 1 + a + In - + 1 - a 3/2 / £ . 50 

1-a + aA 4 ^ (1 - a(l - t)) 2 y/t{t - 1) 

Note that g(X + 0) — g{X — 0) = 27ri on (— (a -1 — 1),0) as this interval is the image of 
the half-axis (— oo, 0) in the z-variable, where it is easy to conclude (cf. (152]) below) that 
g+(z) — g~(z) = 2ni. This jump in the A-variable is also easy to obtain directly from (|50|) . 

Let 

2(1 -af/ 2 f x 1 + «(!-£) dt 



HX) ~ eW2 I { l-a { l-t)) 2 Vt (l-t) (51) 

which is analytic in C \ ((— oo,0) U (l,oo)). Here we choose the branch so that v/i(l — t) 
is analytic in C \ ((— oo, 0) U (1, oo)) and positive for t G (—1, 1). The function h(X) is the 
analytic continuation of g(X + 0) — g(X — 0) off the interval (0, 1). 

Note that 

2 f° 1 + a-2ax , 
h(0) = — — / — dx = 2m. 52 

Now transform the RHP for U as follows: 

f (A) = e -^3/2 C/ ( A ) e -r l (p(A)-i/2) CT3) A G C \ (jO, 1] U {-^^}) • ( 53 ) 

We easily obtain then that T(A) satisfies: 

(a) f (A) is analytic for A G C \ [0, 1]. 

(b) For A G (0, 1) the boundary values of T(A) are related by the jump condition 

T + (A)=T_(A)(^ e Q en , (A) J, AG (0,1). (54) 

(c) T(X) has the following asymptotic behavior as A — > 

f(X) = I + O^X + ^—^y (55) 

Note that the problem is now normalized to I at A = — l=^t. 

Since detT(A) = 1 and T(A) is analytic at infinity, it follows that T(oo) is invertible. 
The function T(A) defined by 

T(A) = f (oo)- l f (A) (56) 
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is the solution to the same Riemann-Hilbert problem as T(A), with the (c) condition replaced 
by 

T(A) — I + 0(1/A), A^oo. (57) 

Clearly, 

f (A) = T- 1 (-*^*\ T(A). (58) 

We now show that the RHP for T is solvable for all < a < 1. For < a < 1 the 
existence of such a T(A) follows simply by pushing forward V(z), the solution of the RHP 
(I35II36I) for the polynomials orthogonal on (0, a) with the weight e~ Anx : the existence of V(z) 
itself follows from the basic results of [201 H2]- So we are reduced to showing that T(A) 
exists in the case a = when the mapping V(z) — > T(A) breaks down. For a = 0, /i(A) = 
41n(vA+z\/l — A), < A < 1. If (A — l) 1 / 2 (resp., A 1 / 2 ) denotes the branch which is analytic 
in C \ [— oo, 1] (resp., C \ [— oo,0]), then in particular (A — 1)+ 2 = —(A — l) 1 / 2 = i\J\ — A, 
< A < 1, and we find 

• (A _l)f + A l/2 - 

(A - I) 1 ! 2 + A 1 / 2 
Thus if r(A) = ((A - l) 1 / 2 + A 1 / 2 )^, then 



2n 



e n,(A) = ^ -r ^ _ (5g) 



e -nfc(A) i \ f(r_/r + ) 2n 1 

e nh W j~ 1 (r + /r_) 2n i ' 



< A < 1. (60) 



Setting Z(A) = T(A)r(A) 2n,J3 , we see that Z(A) solves the RHP: 

(a) Z(X) is analytic for A G C \ [0, 1]. 

(b) For A G (0, 1) the boundary values of Z(X) are related by the jump condition 

Z + (A) = Z_(A)(J \y AG (0,1). (61) 

(c) Z(X) has the following asymptotic behavior as A — > oo: 

Z(A) = (/ + 0(l/A))A nff3 . (62) 



This is the standard RHP for polynomials orthogonal on (0, 1) with the unit weight. There- 
fore the desired solution T(A) exists for a = as well. This completes the proof of solvability 
of the RHP for T(A) for all < a < 1. The above proof of solvability for all n is included 
only for completeness (cf. the last remark at the end of Section 3.5). 

As is standard in applications of the steepest descent method, we now deform the RHP 
as follows. Let S = U^ =1 S J be the oriented contour as in Figure [2j Define a matrix- valued 
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1-cc 
a 



Figure 2: Contour for the S-Riemann-Hilbert problem and the circular neighborhoods £/i i0 
of the points 1, 0. These neighborhoods will be introduced below in connection with the 
construction of parametrices. 

function S(X) on C \ E by the expressions: 

' T(A), for A outside the lens, 



) -^(^) y e -nh{\) \j ' for A in the upper part of the lens, 
T(A) ^ e nh(x) ^ ^ > for A in the lower part of the lens. 

It is easy to verify that S(X) solves the following RHP: 

(a) S(X) is analytic for A G C \ E, where E = Uf =1 E.,-. 

(b) The boundary values of S(X) are related by the jump condition 



(63) 



S + (X) = SI (A) 







AG(S 1 US 3 )\{0,1}, 



where the plus sign in the exponent is on S 3 , and the minus sign, on E 1; 

1 
-1 



(64) 



S + (X) = S-(X) 



X e E 2 = (-1,1). 



(c) 5(A) = J + 0(1/A) as A -> oo. 



For a fixed < e < 1/4, consider the circular neighborhood U\ of radius e at the 
point A = 1. Consider also the neighborhood Uq of A = of radius £3(1 — a) for a fixed 
1/2 > £ 3 > 0. Note that U contracts with growing n for a = 1 — s /(2n) 2 ^. The point 
— (1 — a)/a lies outside U for all a G (0, 1). 
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Figure 3: Form of the contour for the S-Riemann-Hilbert problem. 

In Uq, we can expand the integrand in (}5Tl) in powers of t and t/(l — a): 

2 f x ( 2at ( t 2 \\ 

" W = ' !(0) + ?7VT=^ I (1 + a ~ at) V - — + ° ((T^F ) ) x 

(1 + (/2 + ° (i2) ^^ ,65) 
4v/A / l-6a-3a\ / X 2 \\ 

|A| < £ 3 (1 - a), 

uniformly in a, and where /i(0) = 2iri (see (|52l ). It is the presence of \/l — a in the 
denominator that will allow us to construct a solution to the RHP using a contracting 
neighborhood Uq as a approaches 1. 

We shall now show that the the jump matrices for S(X) on S x U S 3 \ (Ui U U ) are 
uniformly exponentially close to the identity (see (1691) below) as n{l — a) 3 / 2 — > 00. 

To estimate the real part of h(X) outside of the neighborhoods Uq and U\, we now 
describe the form of the lens more precisely. First, we assume that the contour S3 is the 
mirror image of Si, i.e. 

S 3 = Si. 

Therefore, we only need to describe the structure of the contour S x . We assume that for 
< 3?A < 1/2 the contour Sx lies above the straight line originating at zero, and making a 
positive angle 70 with the real axis (see Figure [3]). The value of the angle 70 will be specified 
later on. Similarly, the part of the contour between the vertical line 9ftA = 1/2 and the 
boundary of the neighborhood U\ lies above the line = (1 — 3ftA)tan7! where, again, 
the value of the angle 71 < 70 will be specified later on. Note that the contour S has a 
well-defined limit as a [ 0. 

Let Ao (resp., Ai) be the point of intersection of the contour Si and the boundary of the 
disc Uq (resp., U\) (see again Figure [3]). Let 3?A = /1, SA = u. Thus, A = /i + iu, and on Si, 
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3f£Ao < [i < 3?Ai. Fix some small £2 > 0. Suppose that Si and S 3 are so close to the real 
axis that 

\u\ \u\ 

— < e 2 , < e 2 , KA < a < 3?Ai. (66) 

In particular, this implies that tan 70 < £2 and tan 71 < e 2 . Furthermore, as 1 + a — a/i > 1, 
we have 

\u\ 



1 + a — a/i 

and, as 1 — a + a/i > a/i, 



< w 





< \ u \ 







1 — a + a/i jU 

The above inequalities allow us to perform the following estimate on h(X) for A = fi + iu 
in (Si U S 3 ) \ {E/i U C/ }. Using (EQ), we obtain 

, . , 2(l-a) 3/2 1 + a-au-aiv d(a + iv) 

h(X)=h(ji) + - ' 



e i7r/2 J,, (1 - a + a/i + aiv) 2 ^JJjj, + i v ){\ - ^ - ^) 

,. . 2(1 — a) 3 / 2 1 + a — au /"" / icw \ / zau \ ~ 2 

h(u) + , / 1 H x (67) 

\/fi(l - fi) i 1 - a + OW Jo V 1 + a-a/i/V 1-a + a/i/ 



ic\ 1/2 A if A 1/2 , , , . 2(1 - a) 3/2 1 + a-au r , / i , 

1 + - 1-- dv = h(ji) + -) = - —?-u[l + 0(e 2 )} 

HZ V (1-a + a/i) 2 

where the constant in the error term is uniform for < a < 1. 

The fraction u/(l — a + a/i) 2 in the last equation of fIBTl) can be estimated for some 
£4 > as 

— — > e sin 71 > e 4 , for i < // < 3Ui, 

(1 — a + a/i) 2 2 

|tt| /itan7o tan7o , . 

> 7^ 1 I _.A9 = i h IT7TT\9 > I 68 ) 



^1 — a + a/i) 2 (1 — a + a/i) 2 /i(a + (1 — a)//i) 2 
> £4, for JftAo < /i < -, 



(l + e 2 /(£ 3 sin7o)) 2 " u_^_ 2 , 
where e 4 depends only on e and £j, i = 2, 3, 70, 71, which in turn depend only on e, e 2 , e 3 . 
Since 9Wi(/x) = 0, we obtain from (1681) as n — > 00 for sufficiently small £2 > 0: 

| e -nMA)| = (e-<*), A G Si \ (U U 17 X ), 
| e »fcW| = ( e -Pc) j A G S 3 \ (C/ U U x ) 

uniformly for a G [0, 1 — so/(2n) 2 / 3 ] for some (large) so > and all n > Sq /2 /2, for some 
c = c(e, e 2 , £3) > 0, where 



p = n|l — a 



3/2 
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So except for the jump on the interval (0, 1) and the jumps inside U±, Uq, the jumps of 5(A) 
are indeed exponentially close to the identity as p — > oo. 

For later purposes, we shall need the series expansion of h{\) at A = 0, 1. We have: 

AV\ ( l-6ct-3a 2 / A 2 \\ 

h{\) = 2ixi + — 1 + a + , — A + O -r , A -> 0; 70 

V; eW2^T^V 6(1 -a) V(1-«)V/ 

/i(A) = 4(1 - «) 3/2 v^ (l - (a + l/6)u + (a 2 + 3a/10 + 3/40)w 2 + 0(u 3 )) , A = 1 + u, 
u — > 0. 

(71) 

In ( 170"]) the cut of the root lies to the left of A = 0, and — 7r < arg A < 7r, whereas in (1711) 
the cut lies to the right of A = 1, and < argw < 2ix. 

Note the crucial fact that, as follows from (1701) . ( J7T1) . the quantity n|fo(A)| (resp., n|/i(A) — 
27ri|) is uniformly large on the boundary <9?7i (resp., dUo) for some (large) so > for all 
a G [0, 1 — s /(2ri) 2 / 3 ], if (2n) 2 / 3 > s . Indeed, it is of order s^ 2 for A on dU\ (resp., of 
order n for A on dUo). This will allow us to obtain the desired asymptotic solution of the 
Riemann-Hilbert problem. 

For technical reasons (see the end of the section 4.2. below and also proof of Corollary 
2 in [I7f), we need to control the solution of the RHP for all a G V £o (0) U [e , 1 - s /(2n) 2 / 3 ], 
where T> £0 denotes the disc of radius Eo about zero in the complex a-plane with Eq small. For 
all a G ^0(0) we use ^e fixed contour S = S Q= o in Figure 3 corresponding to a = 0. By 
the preceding calculation we see that \dlh(\; a — 0)| > c > for all A G (S x US 3 ) \ (UqUUx). 
Thus 

| e -nh(A; Q =0)| < e -nc 0; A G (Si U E 3 ) \ (U U Ui). (72) 

Hence, by continuity, we must have 

| e -nfc(A,a)| < g-ncj, (73) 

for all A G (Si U S 3 ) \ {Uq U Ui) and all a G V £o {0), < c' Q < cq, Eq sufficiently small. 

We now begin the construction of parametrices which give, in their respective regions, 
the leading contribution to the asymptotics for the RHP. 



3.2 Parametrix in C \ (U\ U U ) 

First, because of the exponential convergence described above, we expect the following model 
problem to play a role in constructing a parametrix for the solution of the RHP as n — > oo: 

(a) N{\) is analytic for A G C \ [0, 1], 

(b) 

iV + (A)=iV_(A)(_ 1 JV AG (0,1), (74) 
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N(X) = I + O I -J , asA^oo. (75) 

The solution N(X) can be found in the standard way by first transforming N(X) with a 2 x 2 
unitary transformation to the form for which the jump matrix is diagonal and then solving 
the two resulting scalar Riemann-Hilbert problems (cf. [7]). We obtain 

at/as 1 ( m + vnT 1 —i(m — m~ 1 )\ /A— l\ 1//4 

"< A > = 2 U— ■»+«»-' j- m(A) = (— j • (76) 

where m(A) is analytic outside [0, 1] and m(A) — > +1 as A — > oo. Note that det N(X) = 1 
and that N(X) is the unique L p solution of the RHP for any 1 < p < 4. 



3.3 Parametrix at A = 1 

Now let us construct a parametrix in U\. We look for an analytic matrix- valued function 
Pi (A) in Ui which has the same jump relation as S(X) on E fl U\ and instead of a condition 
at infinity satisfies the matching condition on the boundary 

Pi(A)A^(A) = l + 0(l/p), XedU u p = n\l-a\ 3/2 , (77) 

uniformly in A and a as p — > oo. 
Define: 

(e™h(X)/2, for3A>0, 
0(A) " \ h(X)/2, for 3A < • (78) 

This function is analytic in U\ outside (1 — e, 1]. 

We look for Pi (z) in the form: 

Pi (A) = E n (X)P(X)e n ^, (79) 

where E n (X) is analytic and invertible (det E n ^ 0) in a neighborhood of C7x> and therefore 
does not affect the jump and analyticity conditions for P(A)e n< ^ A ) CT3 . 

As Pi (A) is required to satisfy the jump relations (1641) for S, it is easy to verify that 
P(A) = P n (A)- 1 Pi(A)e- n ^ <Ts satisfies jump conditions with constant jump matrices: 

P+(A) = P_(A)(j jY A G ((Ei U E 3 ) n Ui) \ {1}, 

; x (so) 

p + (A) = p_(a)(^_ 1 M, AeE 2 n[/ x . 

Now introduce a mapping of U\ onto a new £-plane 

( = n 2 0(A) 2 = 4n 2 (l-a) 3 M(l-(2«+l/3)M+(3a 2 +14a/15+8/45)M 2 +0(w 3 )), A = 1+u, 

(81) 
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where we used (1711) . The expansion at A = 1 is uniform for a in a bounded set. 

Choosing a sufficiently small e > 0, we see that ((X) is analytic and one-to-one in the 
neighborhood U\. 

Note that if a E [0, 1 — sq/ (2n) 2 ^] then \(\ = 0(p 2 ) uniformly large, if So is large, on the 
boundary dU\ and in a. This is a crucial fact in the present work. When a = 1 — so/(2n) 2 / 3 , 
we have p = s^ 2 /2. 

Let us now choose the exact form of the contours in U\ so that their images under the 
mapping £(A) are straight lines (see Figure H]). Set 

P(A) = *(C), (82) 

So the jump matrices for ^(C) are the same as for -P(A) (they are shown in Figure H]). A 
matrix \I/(C) satisfying these jump conditions was constructed in [27J in terms of Bessel 
functions, namely: 

1) region I 

1 / H^\e~™l 2 e /2 ) H {2 \e-™/ 2 C /2 ) \ 

11/(0 = 2 [nC 1/2 (H^y (e-^l\ l/2 ) < /2 (h®)' (e~^ 2 C 1/2 ) J ' (83) 

2) region II 

1 / H {2 \e™l 2 C 1/2 ) -H^ie^C 1 ' 2 ) \ 

11/(0 = 2 [ -nC 1/2 (tff)' (e i7r / 2 C 1/2 ) vrC 1/2 (i^Y (e^ 2 ^ 2 ) ) ' (§4) 
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*(C) 



5) 



3) region III 

^o(C 1/2 ) ^o(C 1/2 ) \ 

,«C 1/2 ^(C 1/2 ) -C 1/2 K(C 1/2 ) J ' 

where — ir < arg(£) < 7r. 

Here the square root v^C nas the cut on (— oo, 0). Hence, — — n<f)(\) for — n < 
arg(C) < 7r. 

The large-C asymptotics of Bessel functions give (here we choose sq, depending only on 
e, sufficiently large): 



*(0 = ^Vc)- a3/2 



1 z 
i 1 



1 



-1 -2i 
-2i 1 



. -4i 
2 7 C V 4i 1 



uniformly on the boundary 9Z7i. 
Thus 

P 1 (A)= J E n (A)vI/(C(A))e^ A ^, 
where the function E n (X) is found from the matching condition to be 



E n (X) = -j=N(X) 



1 -% 
-i 1 



•3/2 



(86) 



(87) 



Now to complete the construction of the parametrix it only remains to show that E n (X) 
is an analytic function in U\ (clearly, det£' n (A) ^ 0). First, we show that it has no jump 
on the real £-axis. This is easy to verify using the jump condition for N(X) and the identity 
£_ = £+e~ 27rt on the negative half axis. Moreover, a simple calculation shows that E n (X) has 
no pole at A = 1. Thus, E n (X) is analytic in Ui, and the parametrix in JJ\ is given by the 
equations (T79lf82lf83lF84lf85lf88r) for a G [0, 1). 

Below we shall need the first three terms in the matching condition for Pi. Using (186"]) . 
we obtain 

1 



P 1 (A)A^- i (A) = / + Ai(A) + A 2 (A) + O 



X e dU x . 



9) 



Here 
Ai(A) 
A 2 (A) 



8x/C 
3 



N(X) 



N(X) 



2 7 C 



N(X) 



-1 -2i 
-2i 1 



-i 



-i 



-3m 2 + m 2 -i(3m 2 + m 2 ) 
16-v/C V — *(3m 2 + m -2 ) 3m 2 - m" 2 
4i 



2 7 C v~ 4 ^ 



(90) 

where m(A) is defined in (1761) . Note that both Ai(A) and A 2 (A) are meromorphic functions 
in Ui with a simple pole at A = 1. 
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Recall that we use the contour £ = Sq,=o for all a G D eo (0), Eq small. For such a, the 
map A — > C maps Z7i (consisting of the three regions separated by S) onto a set, region to 
region, where the lines separating each region are now no longer straight but lie in small 
cones about the original ones. The opening angles of the cones are proportional to \5sa\. 
Using the same definition for \l/ as in ( 183H851) for each of the new regions I, II, III, we find 
again that fl86l) is valid, and that Pi(X)N~ 1 (X) has the same expansion fl89l) as in the case 
< a < 1. Note that the values of Eq and Sq can be changed (now and below) if necessary. 

3.4 Parametrix at A = 

The construction of the parametrix in Uq is similar. Recall, however, that the radius of Uq is 
£ 3 (1 — a), so it decreases as a — > 1, i.e. as the pole of h(X) approaches the point A = 0. We 
shall see that this neighborhood produces asymptotics for the RHP in inverse powers of n. 

We look for an analytic matrix- valued function Pq{z) in the neighborhood Uq which 
satisfies the same jump conditions as S(X) on £ PI U , and satisfies the matching condition 



P (A)iV- 1 (A) = /+(!-«) 



/2 0(l/n) 



(91) 



uniformly in A on the boundary dU$ as n — > 00. 

Below we define functions in Uq which play the same role as <ft, E n , and P in U\. We 
use the same notation for these quantities as before. Namely, let 




for 3A > 0, 
for SA < ' 



(92) 



This function is analytic in Uq outside [0, £3(1 — a)]. 

As above, we look for the parametrix Pq(X) in the form: 

P (X) = E n (X)P(X)e n ^, 



(93) 



We obtain that 




(94) 



We choose the following mapping of Uq onto a £-plane 



( = n 2 (j>(X) 2 = e~ i7r 4n 2 (l + a) 2 



X 



( 



1 — 6a — 3a 2 



X + O 



1 — a 



1 + 



3(1 - a 2 ) 



where we used ( !70l) . 
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Choosing a sufficiently small £3 > 0, we see that £(A) is analytic and one-to-one in the 
neighborhood Uq. 

Let us also choose the exact form of the contours in Uq so that their images under the 
mapping £(A) are direct lines. In the £-plane the contour and the jump matrices for P(A) 
are the same as in Figure H] with the only difference that all directions are reversed (pointing 
away from ( = 0). It is easily seen that the function 

P(A) = a 3 ^(C)a 3 , (96) 

where \P(C) is given by ( I83H851) satisfies the jump conditions in this case. 
Finally, we calculate E n and obtain 

P (A) = P„(A)a 3 ^(C(A)) ( T3e^ A ) CT3 , (97) 

where 

= ^N(X) Q j)(7rv/Cr 3/2 (98) 

(the analyticity of E n {X) in Uq is verified as above). 
Then we see immediately from (I86|I95I) that 

Pq{X)N-\X) = E^a^iO^^N-'iX) = 1+^0 (-!-] = I + -^=^0 (^j 

(99) 

uniformly in A G 8Uq and a G [0, 1 — so/(2n) 2 / 3 ]. Of course, the bound in fl99|) blows up if 
a — > 1 too rapidly: for < a < 1 — s /(2n) 2//3 , we see that the error term is 0(n~ 2 ^ 3 ). 

Thus the construction of the parametrix in Uq is now complete. 

Using the expansion of ^(C), we can extend ( |99i) to a full asymptotic series in inverse 
powers of n. Substituting ( 1861) into ( |99l) . we obtain in particular: 

P (A)iV- 1 (A) = / + Ai(A) + A 2 (A) + -^=0 (±) , (100) 



where 



1 N(x) f-l 2z\ x _ 1 / m 2 -3m- 2 %{w? + 3m" 2 ) 

l( } " 87C ( j V 2z 1 J iV(A) " TeTJ V^(m 2 + 3m- 2 ) -m 2 + 3m' 2 



(101) 



As above, note that Ai(A) and A 2 (A) are meromorphic functions in Uq with a simple pole 
at A = 0. 

For sufficiently small Eq, the estimate fllOOp extends uniformly for a G T> eo (0) U [0, 1 — 
so/(2n) 2 / 3 ] for all n > as in Section 3.3. 
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3.5 Final transformation of the problem 



Now construction of the parametrices is complete, and we are ready for the last transforma- 
tion of the Riemann-Hilbert problem. Let 



( S(X)N-\X), A 6 C \ (U U Z7i U £), 
R(X) = \ S(X)P \X), Ae{/o\S, 
(S(X)P 1 -\X), A G Z7i \ E. 



(102) 



It is easy to see that this function has jumps only on dU\, dUo, and parts of Si, S3 lying 
outside the neighborhoods Ui, Uq (we denote these parts S out ). The contour is shown in 
Figure [3 Outside this contour, R(X) is analytic. Besides, R(X) = I + 0(1/A) as A — > 00. 

The jumps are as follows: 



R + {X) =R-(X)N(X) 



1 



N(xy 



X e s° ut u s c E 



where the "-" sign in the exponent is taken on S° , and "+", on S 3 
R+(X) = R^(X)P (X)N(X)~ 1 , A G dU \ {intersection points}, 
R + (X) = J R_(A)Pi(A)A^(A)" 1 , A G d\J\ \ {intersection points}. 



out 



(103) 



The jump matrix on S out can be uniformly estimated (both in A and a G [0, 1 — 
so/(2n) 2 / 3 ]) as / + 0(exp(— cp)), where c is a positive constant. In view of the estimates (169]) . 
this is obviously true outside a fixed neighborhood of A = 0, say when |A| > 1/2. However, 
since the parametrix N(X) is of order 1/A 1//4 for A close to zero, and the contour approaches 
A = as a — > 1, we need a more detailed analysis for |A| < 1/2. In that case, we use fl67|) to 
write for all a G [0, 1) (in what follows the same symbols C and c stand for various positive 
constants independent of a, A, and n): 



-n3th(\) 



C 

< —— exp 
v 7 ^ 



-cn- 



vgCj - a) 
1 + afi/ (1 — a))' 



C 



a: 



iexp 



—cn- 



l + at 2 ) 2 _ 



= /(*), 
(104) 
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where t = — a), fi = 3?A. We need to find the maximum value of f(t) in the interval 

*i = a/ (£3/^2) sin 70 < t < / ee t 2 

a/2(1 - a) 

for all a G [0,1 — s /(2n) 2 / 3 ]. For this purpose, it is convenient to consider the following two 
cases separately. 

1) at 2 < 1. Then 1 + at 2 < 2, and we have 

f(t) < ^_ exp[-cnt] = h(t). (105) 
ry 1 — a 



The derivative f[{t) < for t > 0, which implies 



1/3 



max fit) < fxftt) < C^re- cn < Ce7 cn . (106) 



If 1/ y/a > t 2 this is all we need. Otherwise consider 
2) at 2 > 1. Then 1 + at 2 < 2at 2 , and we have 

/(f) < + i— exp[-cn/t 3 ] ee / 2 (f). (107) 
fyl ~~ a 

The only maximum of /2(f) is at the point f c = (3cn) 1//3 . Now choose sufficiently large sq > 
(depending on e 2 , £3). Then 

1 < cn^/s 1 ^ < t, 



Therefore 



max /(f) < f 2 (t 2 ) < Ce~ cn{1 - a)3/2 = Ce~ cp . (108) 



Combining (1106111081) . we finally obtain that the jump matrix on (Si U S3) \ (U U Ui) is the 
identity up to an error of order 



J_ e -nSfe(A) 



Vx 



< Ce- Cp (109) 



for all a G [0, 1 - s /(2n) 2 % 2n > s 3 /2 . 

This estimate can be readily extended to complex a G T> £0 (0). The jump matrices on 
dUo,i admit the uniform expansions given by (1 1 0011891) . 

A consequence of the above considerations is the following result: 
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Lemma 1 Let p = n\l-a\ 3 ^ 2 , a G £> £o (0) U [0, 1 - s /( 2 ™) 2/3 ], U = U \JU lf £ = T, out UdU. 
Also let U\ be the circle centered at A = 1 of radius e/2. Then, for sufficiently small e, 
Ej, j = 0,2,3 (ej, j = 0,2,3 are the e-parameters introduced above in the definition of the 
contour il), there exists s > such that for all a G T> £0 (0) U [0, 1 — s /(2n) 2 ^ 3 }, and n > 
s^ 2 /2, a (unique) solution R(X) of the R-RH problem exists. Moreover, the function R(X) 
admits the following asymptotic expansion, which (and the derivative of which) is uniform 
for a G £> eo (0) U [0, 1 - s /(2n) 2 ^] and all A G U 1} as p — > oo: 

R(X) = 1 + R 1 (X) + R 2 (X) + ■■■ + R k (X) + R { r k+1) (X), 

Ri k+1 \X) = 0{p- k -% ^^ +1) (A) = Oip-x- 1 ), (110) 

k = 1, 2, . . .. The functions Rj{X) = 0(p~ J ) are constructed by induction as follows: 

i?i(A) = ^/ AxOO-^-, R 2 (X) = ^-[ (R l 4s)A l (s) + A 2 (s))^-, (111) 
27T« y ac/ s - A 2m y a{/ s - A 

^(A) = t— / V^-^A,^)^-, i2o=J. (112) 
Z7rz Vac/ ~i s — a 

Remark. The uniformity means that for sufficiently small e, Ej, j = 0, 2, 3, there exist 
positive constants s , c\, and c 2 independent of a, n, A such that 



A /?(*+!) 



< C2 



|#(*+i)| < 

~ ' M ' P fc+1 (113) 

VA G f/i, Va G P eo (0) U [0, 1 - s /(2n) 2/3 ], Vn > Sq /2 /2. 

We also note that, p > s 3 /2 /2, Vet G £> eo (0) U [0, 1 - s /{2n) 2 / 3 ], n > s 3 Q /2 /2. 

Proof of Lemma 1. We shall follow a similar line of argument to the one which was 
used to prove similar statement in [17] (Lemma 1). For simplicity, as in [T7J, we will only 
prove the expansion (IllOp in the case k = 2, which is all that is needed for the problem at 
hand. We shall also adopt the notation: 

R^(X)=R r (X). 

Besides, as before, the symbol c will stand for various positive constants independent of a, 
A, and n. 

Write the jump condition for R(X) in the form 

R + + Ri+ + R 2 + + R r + = (Ro- +Ri- + R 2 -+ R r -)(I + A 1 + A 2 + A r ). (114) 

Here Aj and A 2 are given by fll01ll90p on OUq, dU\, respectively, and we set Ai = A 2 = 
on the rest of the contour. A direct analysis of the expressions fll01|90l) shows that A k = 
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0{{n- k \l - a]- 1 / 2 ) on dU , and A fc = 0{p- k ) on dU^ Similarly, A r = 0(1/ p 3 ) on dU x (this 
error term arises from the Bessel asymptotics) , A r = 0(|1 — a| 4 /p 3 ) on dU , and, by (1 1 09j) . 
A r = 0(e~ cp ) on E \ <9C7. 

We now show that we can define i?i and R2 so that they are of order 1/p and 1/p 2 , 
respectively. We then show that the remainder R r is of order 1/p 3 . Set 

Ro = I- 

We define Rj by collecting in (11141) the terms that we want to be of the same order. First, 

R 1+ (\) = R 1 -(\)+A 1 (\), AgE. (115) 

We are looking for a function i?i(A), which is analytic outside S, satisfying i?i(A) = 0(1/A), 
A — > 00, and the above jump condition. The solution to this RH-problem is given by the 
Sokhotsky-Plemelj formula, 

Ri(X) = C(Ai), (116) 

where 



2ni Jf, s — A 

is the Cauchy operator on E. The condition Ai(A) = 0(1/ p), A G S, p — > 00 (uniform in 
a), implies that there exist c, 5, so > such that 



3/2 

\Ri{X)\<c/p, n>^- (117) 



uniformly in a G £> eo (0) U [0, 1 - s /(2n) 2/3 ] and A satisfying dist(A, E) > 5. Actually, this 
estimate is uniform for all A G C \ E up to E. Indeed, since 

RlW = 7^-1 Ms) -^r, (118) 
2m J du s - A 

for A outside a fixed neighborhood of zero, this is seen by shifting the contour to a fixed 
distance from the point A. Inside that neighborhood, the distance of the shift will depend 
on a. Namely, the distance is e'\l — a\ for a fixed (sufficiently small) e' > 0. Then 

|C(A0| < max lA^— — + - < ° +- = -, (119) 
c|l — a\ p n\l — a\ 6 ' z p p 

on and close to OUq. Here we used the estimate Ai = 0(n~ 1 \~ 1 / 2 ), so that in the neighbor- 
hood of the circle dUo the inequality 

c 

max I Ai| < — 
n v 1 — a 
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holds. It should be observed that, by the same deformation of the contour of integration in 
(11181) . one obtains the analytic continuations of both the functions i?i + (A) and i?i_(A) in the 
neighborhood of the contour dU and hence in the neighborhood of £ (we note that on the 
part S out of the contour £ R\(X) has no jump). Moreover, the estimate (11 17ft is preserved 
under this analytic continuation. 

Now define -R 2 (A) by the jump condition 

i? 2 +(A) =i?2-(A) + i?i_(A)A 1 (A)+A 2 (A), A G S, (120) 

together with the requirement of analyticity for A G C\S, and the condition i? 2 (A) = 0(1/A) 
for A — > oo. The solution to this RHP is 

R 2 (X) = C(i? 1 _A 1 + A 2 ), AgC\£. (121) 

Using ( I117p and the estimates for A 2 , we obtain in the same way as for Rx, 

3/2 

|^ 2 (A)| <c/p 2 , AgC\£, n >\ ( 122 ) 

with the same uniformity and analyticity properties in a and A. 
Now from (1114011501201) we obtain 

R r+ (\) = R r _(\) +M(A) +# r _(A)A(A), A G E, (123) 

where 

M = R 2 -^i + {Ri~ + i? 2 -)A 2 + (/ + #1- + R 2 -)\, A = Ai + A 2 + A r . 

Remark. In the terminology of p2], equation ( 1 1 2 3 f) is an inhomogeneous RH-problem 
of type 2. 

Since R r = R — I — Rx — R 2 , the matrix function i? r (A) is analytic outside S and satisfies 
the condition R r (X) = 0(1/ X) as A — > 00. Therefore 

i2r(A) = C(M) +C( J R r _A), AgC\S. (124) 

Hence 

R r _(\) = C-(M) +C_(i? r _A), AGS, (125) 

where C_(/) = liniA'^A C(/), as A' approaches a point A G S from the — side of E. Now 
defining the operator 

C A (f) = C-(fA), 

we represent ( 11251) in the form 

(I-C A )(Rr-)=C-(M). (126) 
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By virtue of the estimates (11001) . (!89|) . and (I109P we have that 



l|A|lL 2 (S)nL°°(s) ^ ~i (127) 

for all a G V eo (0) U [0, 1 - s /(2n) 2 / 3 ] and n > s 3 J 2 /2. 

The Cauchy operator C_ is bounded in the space £ 2 (X) (see, e.g., |28j), and by a 
standard scaling argument (the Cauchy operator is homogeneous of degree 0), its norm is 
bounded by a constant independent of a. This together with the L°° part of the estimate 
(11271) implies that the operator norm ||Ca||l2 = 0(l/p), and hence I — is invertible by 
a Neumann series for sq (and, therefore, p) sufficiently large. Thus (I126p gives 

R r - = (I-C A )-\C4M)), (128) 

and this proves the solvability of the -R-RH problem for all a G T> £0 (0) U [0, 1 - s /(2n) 2 / 3 ] 
and n > Sq/2. Moreover, using the L 2 part of the estimate (|127p . we conclude that 
||C_(M)|| L 2(2) = 0(p~ 3 )- Together with (I128P this yields the uniform estimate 

||-Rr-|| L 2(E) < ^3, (129) 

Va E V £o (0) U [0, 1 - s /{2nf/% n > s 3 Q /2 /2. 
The solution -R(A) of the i?-RH problem is given by the integral representation 

R(X) = I + R 1 (X) + R 2 (X) + C(M) + C( J R r _A)(A), (130) 

A G C\S. 

Remark. Let Q^, k = 1, 2, 3, 4 denote the connected components of the set C\S. Then, 
using again the possibility of the contour deformation when solving the integral equation 
), and taking into account the triviality of the jump matrix monodromy at each node 



point of the contour S, we conclude that the restriction i?|n fc (A) is continuous in for each 
k (see e.g. [1]). This means that equation (11301) defines the solution of the i?-RH problem 
in the classical, point-wise continuous, sense. 

Combining the inequality (11291) with equation (I130p . we can complete the proof of the 
lemma. Indeed, assuming that A G U\, we immediately obtain the estimate 



\C(M)(\)\<-, n>sf/2, (131) 



for the fourth term in the r.h.s. of (11301) . and the estimate 



|C(i? r -A)(A)| < c\\R r .\\ L2{t) m\ LHt) < iL (132) 
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n > sf/2, 

for the fifth term. Both the estimates are uniform in a G £> eo (0)U[0, l-s /( 2 ™) 2/3 ]- Together 
they yield the estimate 



c 



|i2r(A)|<-^, n>sl /2 /2, (133) 

uniformly in a G D eo (0) U [0, 1 — so/(2n) 2//3 ] and A lying in Ux- This establishes part of the 
estimate (IllOp for the error term. The estimate for the derivative follows immediately from 
( 11241) . This completes the proof of the lemma (in the case k = 2). □ 



Remark. (Cf. Remark 2 in [TT].) Part of the assertion of Lemma 1 is that the solution 
of the i?-RH problem, and hence of the original T-RH problem, exists and is unique for all 
a G T> eo (0) U [0, 1 — s /(2n) 2 / 3 ] and n > s^ 2 /2 with s sufficiently large. This is all we need 
in the analysis that follows; however, the solution of the i?-RH problem actually exists and is 
unique for all a G T> £0 (0) U [0, 1) and alln > for some (possibly smaller) e > 0. Indeed, by 
the discussion following ([58]) . the T-RH problem, and hence the i?-RH problem, is solvable 
for all a G [0, 1), n > 0. Since, by the previous remark, the solution of the i?-RH problem 
is continuous up to the contour, the problem is easily seen to be solvable for a G V e / (0), 

< n < s^ 2 /2 for some e' Q > by continuity of the jump matrix at a = 0. By Lemma 1, 
the i?-RH problem is solvable for all a G V £o (0), n > s^ 2 /2. Thus the i?-RH problem, and 
hence the T-RH problem, is solvable for all n > on P e "(0) U [0, 1), where e" = mm{e , e' }. 



4 Evaluation of the differential identity 



4.1 Exact transformations 



We start with the differential identity (I3TI) . Note that since V(z) is related to Z7(A) by the 
expression (I46|45l) 

17(A) = V(z(X)), z 



aX 



1 — a + a\' 



we have 



dX 

dz 



all — a) 



and (|37|) can be rewritten in terms of Z7(A) as follows 



4- InDJa) 



-ina 



da " 2ma(l — a 

Note that the derivatives in ( 11351) are taken w.r.t. A 



(U U (1)U^(1)-U' U (1)U 21 (1)). 



(134) 



(135) 
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By (I53|58p . the matrix elements of U(X) can be expressed in terms of T(X) as follows: 

1 — a" 



U U (X) 
U 21 (X) 



T -l 



a 
1 — a 



a 



T(A) 
T(A) 



o n 9( x ) 



11 



(136) 



e -nl e ng{\) 



21 



Furthermore, for A outside the lens in U\ 

T(X) = S(X), S(X) = R(X)P 1 {\). 

Note also that by (JS7D 

= ( J R(A)E n (A)vI/(C)) J ie^ A ), j = 1, 2, 
and, as follows from the definitions of the functions 0, ft,, and the properties of g(z), 

2z(X)+l/2, 



(137) 



0(A) + = + 5 = =F ^ + 0± = — 



where g±(A) stand for the analytic continuation of these functions. Here the upper sign 
corresponds to SA > 0, and the lower, to 3 A < 0. 

Hence, (I135[) finally gives 

(1 ' ' ' ■((iiE*) 11 (l)(iiB*) / 21 (l)-( J R£;*) , 11 (l)(iiE*) 21 (l)), (138) 



da ^ 2iria(l — a 

where we used the fact that detT _1 (— (1 — a) /a) = 1. In f 1138D . the derivative at A = 1 is 
taken along a path in U\ outside the lens. In the next subsection we use the solution of the 
Riemann-Hilbert problem for R(X) (found in Section 3) to construct the asymptotics of the 
r.h.s. of ffT38|) . 



4.2 Asymptotics 



Consecutive asymptotic terms in the expansion of the logarithmic derivative (11381) are gen- 
erated by consecutive terms in ( 11101) : 

R(X) = I + R 1 (X) + R 2 (X) + ■ ■ ■ . 

Thus, setting R = I in (11381) gives the main asymptotic term of 4- In D n {pt): 



1 



-((^)n(l)(^) 21 (l) - (^(^^(l)), 



2iria(l — a) 
Using (1581) and (185|) . we obtain 

(^)n(C) = MA), (S*) 2 i(C) = -^-(A), 



(139) 



(140) 
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where 



A*±(A) = \I^C 1/4 (rn-\X)I (^C)±m(X)Il ) ( v ^)). 



(141) 



Using the expansion of Bessel functions as ( — > (i.e. A — > 1), we obtain 
H±(l) = M = V5m(l -«) 3/4 , /4(l) = a±6, 



a = M 



n 2 (l-af 



a 



6 



Mn(l - a) 



3/2 



(142) 



Substituting these values into (11391) . we find the main asymptotic term 



d 

da 



In DJoc 



n \ 
~ — (1 — a) 

a 



To obtain the next term, we need to compute first 



fli(l) 



2m 



du 



Ax (A) 
A - 1 



it^l) 



2vu 



<?(/ 



Ax (A) 
(A -1)2 



tZA. 



(143) 



(144) 



We now examine A 12 in the neighborhoods of the points A = 0, 1. Using (195]) and expanding 
the matrix elements of N(X), we obtain from (11011) : 



F (a) 



F (a)-5/2 z(F (a) + 7/2)\ 

A ^ 32n(l + a) V*(Fo(a) + 7/2) -F (a) + 5/2 ^ ^ ^ {Ah 

l-6a-3a 2 „ v 7 !"^ /-I 



32n(l + a) \ -i 1 



(145) 



Ae % 



1 4z 
-4z 1 



Aef/„ 



(146) 



6(1 -a 2 ) ' 
For A 2 (A), we obtain similarly: 

= 3(1 -a)(l+ 0(A)) 
21 ' 2 9 n 2 (l + a) 2 A 

In Ui, a similar calculation based on (190|) and (IHTl) gives (A = 1 + it) 

A m = — -I - ( - 5 /2 + a + l/6 -i(7/2 + a + 1/6) \ 

ll ' u 32n(l -a) 3 / 2 V-^(7/2 + « + l/6) -(-5/2 + a + 1/6) J 

u f 3/2- (5/2) (a + 1/6) + F x (a) -i(-3/2 + (7/2) (a + 1/6) + F x {a)) 

32n(l - a) 3 / 2 V -i(-3/2 + (7/2)(a + 1/6) + Fx (a)) -(3/2 - (5/2) (a + 1/6) + F x {a)) 



+ 



+ 



0(n 2 ), 
FAa) 



a 



1 + u = A, 
31 



6 4 • 45 ' 



Aet/i, 
A 2 (A) = 



1 -i 
-i -1 



32n(l - a) 3 / 2 
3(1 + (2a + l/3)w + 0(u 2 )) ( -1 4i 



2 9 n 2 (l -afu 



-Ai 



(147) 
A g fA. 
(148) 
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Now the expressions for R(l) and R'(l) are obtained from the above results and (11441) 
by a straightforward residue calculation: 



5 7] 

r\ —5 



32n(l - afl 2 

a t 
t —a 

% 



7/2 + a + 1/6 



32n(l-a) 3 / 2 _ 
(1-a) 



5/2 -or -1/6 - 



1 + a 



1 + a 



(7 



32ra(l - a) 3 / 2 



32n(l - a) 3 / 2 
-3/2+(7/2)(a+l/6) + Fi(a) + 



3/2 + (5/2) (a + 1/6) - Fx (a) 
(l-«) 2 " 



1 + a 



1 + a 



(149) 

Note that the contours dUo t i are traversed in the negative direction. 

We shall be using the following notation for the expansion terms of the logarithmic 
derivative ( 1138p . We denote Rk ■ R m (Ro = I) the term given by 

— — \ ((it! fc ^) 11 (l)(it! m ^) / 21 (l) + (R m E^) n (l)(R k E^)' 21 (l) 

2ma{l - a)(l + d k>m ) (150) 

-(R k E^)' n (l)(R m E^) 21 (l) - (R m E^)' n (l)(R k E^) 21 (l)). 



For example, the main term (I139P is I ■ I. We can now evaluate the next (Ri ■ I) term in the 
expansion. It is written as follows: 



— hi £„(«)- -(1-a) 2 --—- A\Rx{\) 1 

da a 2ma{l — a) \ { 

M ( -%% ) }[ w - ( -%% ) }', 



M(-i//_(l))+ 



(151) 



^(1) ^(1) 



—i 



M 



a 



M 2 (r + ia) 
iria(l — a) 4(1 — a 2 ) 



where we first simplified the expression substituting the above symbolic representation of Ri 
in terms of 5, r], a, r, and used their numerical values only at the last step. 

It turns out that the two terms in the asymptotics just obtained is all we need (up to 
the error term). The following lemma is the main result of this section: 



Lemma 2 There exists sq > such that the expansion 

d TL 2 Oi 

— In D n (a) = — {l-af + — — + r(n,a), 

da a 4(1 — a 2 ) 



r(n, a) 



rM 1 

1 — a \p 



p = n 



l-a| 3/2 , 



(152) 
(153) 
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holds uniformly in a 6 (0, 1 — sq/ {2n) 2 ^\ for all n > s^ 2 /2. 

Proof. It only remains to prove the expression for the error term. We consider the 
expansion of R(X) up to the third term: R = I + R\ + R 2 + R r - Since Rk = 0(p~ k ) and, 
according to ffT4"2j) . /4(1)//±(1) = 0(p 3 ), it is not difficult to deduce from ffT38|) (cf. (jI5Il) ) 
that the contribution of the terms R r ■ Ri, R 2 - R 2 and higher are of order (a(l — a))~ 1 0(p~ 1 ). 
Thus we shall need to consider in detail only the following 4 terms: R\- R±, R 2 ■ I, R% ■ Ri, 

Rr-I. 

For the Ri -R\ term, which we denote Lu, we obtain after a calculation similar to fll 51 j) : 

n 2 1 

L n = (1 - a) 2 {5 2 + rj 2 ) = — — -(a + 2/3)(2 + 5a - a 2 ). (154) 

a 2 8 a(l — a^J 

For further analysis, we need to calculate i?2(l)- It is given by the formula: 

= 1 / fll -(A)A (A) + A 2 (A) ^ 

2m y 8i7 a - 1 

The solution of the Riemann-Hilbert problem for Ri inside U\<a is given by the expression 
(which we write on the boundary) 

i?i-(A) = ^ i T + y-A 1 (A), A G dU, (156) 

where Ax, C\ are defined in fll47|145l) . Note that outside C/^o the solution is 

It is easily seen that the jump, analyticity conditions, and the condition at infinity of the 
Riemann-Hilbert problem for Ri(X) are satisfied, and therefore, by uniqueness, this is the 
solution. 

The expansions for A 12 obtained above and the formulas (I156|ll55p give, by a residue 
calculation, the final expression for i?2(l): 



/3 7 J ' 2 9 n 2 (l-a) 3 



1 3(1 + a)/ 1 + a 



57) 

where the expression for (3 is omitted as it is not needed below. 

To compute the " R% • I" term (which we denote L 2 o) note first that the contribution of 
the terms in that expression involving -R' 2 (l) is of order (a(l — a)) _1 0(p~ 1 ) and we need not 
calculate them. The remainder gives a nontrivial contribution, and we obtain: 

L^^- )+ -J— Of-T*). (158) 
7ra(l — a) a[l — a) \n{l — ay i '' ! ' J 
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The expression for 7 tells us that this is equal to —L\\ (11541) up to the error term. Thus, we 
conclude that the contributions of R2 ■ I and Ri ■ R\ terms cancel each other. 

The analysis of the R2 ■ R\ term is now easy to carry out, and we find that this term is 
of order (a(l - a))~ l O(p- 2 ). 

For any matrix elements of R r (l) (we only know they are of order 0(p~ 3 )), we obtain 
that the R r ■ I term is of order (a(l — a))~ l O(p~ l ). 

Thus, in view of uniformity of the error term in the expansion of R(X), the lemma is 
proven but with the remainder 

r(n, a) = — - — -O ( - J . (159) 



a(l — a) \p 

We now show that a in the denominator here can be omitted. First, we notice that r(n, a) = 
O n (l) as a — > and n is fixed: this follows immediately after substitution of the expansion 
(1271) into the l.h.s. of (I152p . However, we need an estimate which is uniform in n. To obtain 
such an estimate, we use the extensions of our expressions for complex a discussed above. 
As follows from (115211271) . r(n, a) is an analytic function of a in T> eo (0). Thus 

r(n,a) = — f ^-^d5, Id < e„/4. (160) 

2m JoD eo/i {o) « - « 

Since by (11591) . r(n, a) is uniformly bounded on &D eo / 2 (0), it follows that r(n, a) is uniformly 
bounded by 0(1/ p) for all a e P £o / 4 (0), and all n > s^ 2 /2. Lemma 2 is proven. □ 



5 Proof of Theorem 1 

Integrating the differential identity (1 152f) from a (close to zero from above) to any «o < 
a < 1 — So/(2n) 2 / 3 , we obtain: 

ln J D n (a)-ln J D n (a ) = n 2 In 2 (a - a ) + — — --In- ^ + 



«o V ' 2 ; 8 1 - a 2 Q \n(l - afl 2 

(161) 

for all n > s^ 2 /2. Note from (11531) that the term 0(l/n(l — a) 3 / 2 ) does not depend on a - 
Substituting for lnD n («o) the expansion (j27jl and taking the limit «o — ^ 0, we obtain for 
any < a < 1 - s /(2n) 2 / 3 , 



In D n (a) — n I — h In a — 2cH J Inn ln(l — a )+ 

i / i V 12 8 (162) 

ln2 + C'(-l) + 0(— r^)+<y„. 



12 s v ' \n{l-a) 3 / 2 
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Fix any s > So and, for n sufficiently large, set a = 1 — s/{2n) 2 ^. Now take the limit 
n — > oo. As n — > oo, the r.h.s. of f 1 1 6 2 j) becomes 

- ^ - I In 8 + ± In 2 + C'(-l) + 0( S - 3 / 2 ). (163) 

On the other hand, as s is any fixed number s > sq, the l.h.s. of (11621) converges to 
lndetCT-iQ by pT]). □ 
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6 Appendix 

Here we present an alternative derivation of the identity (|37|) . Let 

(f>(x) = -uj n -i(x), i/;(x) = ~uj n (x). (164) 
The determinant ffTTj) is written then as follows: 

D n (a) = det / X(a,oo) ■ (165) 

The operator K(x,y) = (0(x)V'(l/) — (p(y)ip(x)) / (x — y) is of integrable type, and hence 
(see, e.g., [23l El [9]) D n (a) is related to the following Riemann-Hilbert problem for a 2 x 2 
matrix- valued function Y(z) (Figure [6]): 

(a) Y(z) is analytic for zgC \ [ot, oo). 

(b) Let x G (a, oo). Y(z) has L 2 boundary values Y + (x) as z approaches x from above, 
and Y_(x), from below. They are related by the jump condition 

v t \ v ( \ t \ i \ f 1 + 27ri0(x)V'(a;) -2iii^(x) 2 \ , 

Y + (X) = Y^X)V Y (X), V Y (X) = ^ l _ ), X 6 («, 

(166) 
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+ 

a 1 

Figure 6: Contour for the y-Riemann-Hilbert problem, 
(c) Y(z) has the following asymptotic behavior at infinity: 

Y(z)=I + oO^\ , asz-^oo. (167) 

As in [SI [21], it is possible to reduce the RHP for Y to an equivalent RHP with an 
"elementary", in fact constant, jump matrix (see (11740 below). Note first that for any 
functions ip(x), cj)(x) such that ij;(x)(f)(x) — <f>(x)ip(x) = 1, we have 

M.) = A(*)(J "f )A-H*), *<*>), (168) 

Note that the condition on if)(x), <j>{x) is equivalent to the following one: 

detA(x) = 1. 

Let 

2nz f°° </>(€) c -2n 



^(z) e 2 " 2 / °° |§e- 2 <^ 



*(*) = K. ^rocfcl^,,; • (169) 



The function $(z) is analytic in C\R+. Using the orthogonality property of the polynomials 
p n (x), p n -i(x) with respect to the weight e~ Anx , we see that $(2) solves the following RHP 
on E + : 

(a) $(2) is analytic for zgC \ [0, 00). 

(b) For x G (0,oo) the L2 boundary values $+(x) and $_(x) are related by the jump 
condition 

*+(*) = S-(*)fi 2 fV xG(0,oo). (170) 



x 1 

(c) $(2) has the following asymptotic behavior as z — > 00: 



*(:): (/ + ^(;))(f^ n ) ff3 - (171) 
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By standard arguments, see [7], det $(z) = 1. Hence, we see that for x > 0, we can take 
in (PHD 

A(x) = (172) 

The decomposition (11681) suggests the following transformation of the Riemann-Hilbert 
problem. Let 

X(z) = Y(zMz), (173) 
It is easy to verify that X(z) satisfies the following problem: 

(a) X(z) is analytic for z G C \ [0, a}. 

(b) For x G (0, a) the L 2 boundary values X + (x) and X_(x) are related by the jump 
condition 

X + (x) = X_(x) 2 ™), xG(0,a). (174) 

(c) X(z) has the following asymptotic behavior as z — > oo: 

= (l + O (iY\ f^ e -^V 3 , (175) 



Thus satisfies the same RHP as 3>(z), but now on the interval (0, a). 
The transformation 

V(z) = (v^y) X (^) e2nZ<T3 (^T' s/2 (176) 

converts the RHP to the RHP for V(z) of Section 3. 

We now turn to the derivation of the identity for D n (a). Write the determinant (11651) 
in the form 

D n (a) = det(I - K), 
where K is an integral operator acting on functions f(x) from L 2 (a, oo) as follows: 

(Kf)(x)= K(x,y)f(y)dy, K(x,y) = 



x-y 

The logarithmic derivative of D n (a) w.r.t. a has the form 

^A,(a) = - t r ((I-K)-^) = ( ( ,-K^K )(a , a)= m 
((/ - K)-\K -1 + I))(a, a) = R(a } a), 
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where R(x, y) is the kernel of the operator (J — K) 1 — I. As noted above, the kernel K(x, y) 
has the structure of an "integrable" kernel. A consequence of this fact is the identity 

R{ ) = z FM3M±3(*mi, (m) 

x-y 

where the Fj(z) are expressed in terms of the solution of the Riemann-Hilbert problem for 
Y(z) as follows: 

Fj{z) = Y +J1 iP + y +ii2 0, j = 1, 2. (179) 
Comparing this with the definition (11731) of X(z) we see that 

F J (z)=X jl (z), j = 1,2. (180) 

Substituting then R(a, a) = ]im. x _> a R(x, a) into (j 1 77j) . we obtain: 

^\nD n {a)=X u (a)X^(a)-X' n {a)X 21 (a), (181) 

which expresses the logarithmic derivative of D n (a) in terms of the solution of the Riemann- 
Hilbert problem for X(z). Now the function X(z) is related to V(z) by the expression (j!76p . 
In particular, 

X u (z) = ^e-^Vniz), X 21 (z) = -^— e - 2nz V 21 (z). 

Calculating the derivatives of these quantities at z = a and substituting into (11811) . we finally 
obtain (l3j 
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